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' Abstract 

I We summarize the theory and phenomenology of the Color Glass 

Condensate reviewed previously by E. lancu and the author in hep-ph/0 3032C)4| 
Qh! In addition, we discuss some of the subsequent developments in the 

^ • past year both in theory and in phcnomcnological applications. 



■ 1 Introduction 

In the spring of this year, I delivered three 90 minute lectures on the Color 
Glass Condensate (CGC) at the joint Hadron 2004/RANP meeting of the 
Brazilian high energy and nuclear physics communities held at Angra dos 
Reis in the state of Rio De Janeiro, Brazil ^. A sizeable fraction of the 
material covered there was discussed previously by E. lancu and myself in 
a review last year . I will therefore only summarize those topics and refer 
the reader to Ref. for details. This summary may be useful to the reader 
who might be interested in a quick review of the subject. Subsequent to the 



'^Subsequently, a version of these lectures was presented at HUGS 2004 at Jefferson Lab, 
and at the University of Bielefeld. Part of the material discussed here was also presented 
at the CTEQ summer school in Madison, Wisconsin and at the Erice conference QCD at 
cosmic energies. 
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review in Ref. PP, there have been several important developments in both 
the theory and phenomenology of the CGC. (The latter in particular has 
been greatly stimulated by the remarkable data on dA scattering taken by 
the RHIC experiments [2|.) Some of these developments were discussed at 
Angra by myself as well as by lancu |H] and McLerran |^. We shall discuss 
these developments wherever appropriate. Hopefully these lectures will then 
present an up to date perspective on the CGC for interested readers. A few 
other reviews which overlap with the material covered in these lectures are 
listed in Ref. [H]. 

1.1 Summary of Three Lectures on the CGC 

The outline of these lectures is as follows. In Lecture I, we begin with a gen- 
eral introduction to outstanding problems in high energy QCD. We intro- 
duce the paradigm of Deeply Inelastic Scattering (DIS) and in this context 
discuss the evolution equations of perturbative QCD. A likely consequence 
of QCD evolution at small x is saturation of the wavefunction [^1 [7j in the 
limit ^ of = fixed and s ^ oo. This leads to the formation of a Color 
Glass Condensate [HIE, for < Qs(x), where Qs is the saturation scale. 
We discuss a simple model of DIS, the Golec-Biernat-Wusthoff model IQ], 
which provides a saturation inspired phenomenology of DIS. 

In Lecture II, we begin with a discussion of key features of Light Cone 
quantized QCD. In particular, the parton model picture of QCD is manifest 
in Light Cone quantization. It provides a consistent framework to discuss 
the multi-gluon components of the infinite momentum frame wavefunction 
that dominate high energy scattering in QCD. This discussion provides the 
background and the motivation to construct a classical effective field theory 
that describes the ground state properties of hadrons at high energies [5]. 
We next discuss the small x renormalization group equations (often called 
JIMWLK equations) which consistently include quantum corrections to the 
classical effective theory as one goes to higher and higher energies The 
full structure of solutions to the JIMWLK equations is extremely compli- 
cated and only mean field analytical solutions are known 11 . A simplified 
equation which captures some of the non-linear structure of the full RG 

^Here = — > 0, where is the momentum transfer squared in deeply inelastic 
scattering, and s — [P + fc)^ is the center of mass energy squared in the collision of an 
electron with four momentum k and a hadron with four momentum P. The Bjorken 
variable x is defined as a; = Q^/2P ■ q, and one has xy ^ Q /s, where y = P ■ q/P ■ k 
is called the inelasticity. For simplicity, we will assume throughout that y = 1 and hence 
X « Q^/s. 
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evolution is the Balitsky-Kovchegov (BK) equation The BK equation 
has not been solved analytically to date though numerical solutions of these 
equations exist ^1 EI • These are discussed briefly. The BK equation 
can be solved in a "diffusion" approximation ^in analogy with the theory of 
travelling wave fronts in statistical physics j64j . The analogy is a powerful 
one and very recently it has led to interesting developments in the treatment 
of fluctuations at high energies 1^ IBS] • 

Lecture III deals with the application of the CGC picture to describe pp, 
p/D-A collisions and A-A collisions. We begin by briefly describing theoret- 
ical attempts in the CGC /saturation picture to fit the HERA DIS data. We 
next briefly discuss the validity of kj_ factorization in hadronic scattering at 
high energies. This is of phenomenological interest because some saturation 
models employ k± factorization and others do not. The recent RHIC data 
on Deuteron-Gold collisions has remarkable features especially at forward 
rapidities. We discuss the CGC interpretation of this data and argue that 
the unusual features of the data at forward rapidities have a natural expla- 
nation in the CGC framework. Finally, we discuss CGC based approaches 
to the phenomenology of heavy ion collisions. While there is an emerging 
consensus that final state interactions are essential to understand the 
RHIC data, the initial state (as described by the CGC) helps constrain the 
properties of the bulk matter formed in heavy ion collisions. An outstanding 
theoretical problem is whether thermalization is achieved starting from the 
CGC. There has been much recent excitement about this possibility and we 
shall briefly summarize that discussion. 



2 Lecture I 

Perturbative QCD is very successful but it describes only a very small part 
of the total cross-section. Recall that 

^Rutherford g2 ' ^'^^ 

which is negligibly small in the limit of oo. The bulk of the semi- 

hard and soft cross-section remains to be understood. Lattice QCD is a 
first principles approach but it has been successfully applied so far only to 

■^It was later argued by these authors (in the second and third papers of Ref. 64 ) 
that the analogy with travelling waves goes beyond this approximation and is generic of 
saturation models in QCD. This is because in the latter, like the former, is characterized by 
unstable, rapid linear growth that is damped by non-linear effects which lead to saturation. 
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static quantities. Its applicability is least at high energies where higher mo- 
ments of leading twist operators become as important as lower moments, 
and moreover, where higher twist operators begin to contribute significantly 
to physical processes. The fact that the two most rigorous techniques, per- 
turbation theory and the lattice, are inadequate is especially unfortunate 
because it appears that there may exist simple structures controlling the 
high energy behavior of cross-sections. These regularities have been known 
for a long time and understanding these in the language of the fundamental 
theory remains an outstanding scientific puzzle. One such example is the 
parametrization of total cross-sections as 



advocated by Donnachie and Landshoff [20] . The first term corresponding to 
rising cross-sections with energy is believed to occur due to the exchange of 
a particle with vacuum quantum numbers (the Pomeron), while the second 
term (corresponding to falling cross-sections) is due to Reggeon exchange. 
Such simple fits work extremely well for a large number of observables. It 
must be noted that alternative parametric forms are argued to work just as 
well. Assuming that the Donnachie-Landshoff fits, by the logic of Occam's 
razor, correspond to reality, what are Pomerons and Reggeons and how can 
we construct them from the underlying theory? The most sophisticated at- 
tempts to understand the Pomeron involves the exchange of two "reggeized 
gluons" -the Pomeron thereby constructed is called the BFKL Pomeron j21j . 
The BFKL Pomeron is constructed in perturbation theory (very sophisti- 
cated perturbation theory) but it has a few problems. More sophisticated 
work can cure some of these problems but the fact remains that the BFKL 
Pomeron does not describe total cross-sections and indeed is hard pressed 
to explain perturbative results as well. 

A "road map" of the strong interactions shown in Fig. ^ is useful at this 
stage. On the y-axis, we plot the Rapidity (or the Energy) and on the x-axis 
the momentum resolution (Q^). At low and low energies, Chiral Pertur- 
bation Theory works very well. Meson exchange modes do a reasonable job 
at somewhat higher energies (but still small values of Q^). At high energies, 
and small momentum transfers, Pomeron phenomenology is successful if not, 
as discussed, understood. On the other extreme of the plot, the perturbative 
QCD evolution equations work very well. In perturbative QCD, one obtains 
large logarithms in asln{Q'^ /QQ)ln{xo/x), where xq and correspond to 
the initial values in the evolutions. The DGLAP equations ^22^ compute the 
leading logarithms in ln{Q'^/QQ), while the logarithms in ln(xo/a;) are 
sub-leading. The situation is reversed for the BFKL equation, where the 



4 




Momenlum Resolulion Q'3 



Figure 1: A road map of the strong interactions. 

leading logarithms are those in asln(xo/x). Both approaches lead to lin- 
ear evolution equations, in ln((5^/(5o) and ln(xo/x) respectively, and both 
lead as well to rapidly rising cross-sections at small x. Very rapidly rising 
cross-sections violate unitarity and there must therefore exist a mechanism 
in the theory which restores unitarity at small x. Since the unitary bound- 
ary may be reached even for » Aqqj-, at sufficiently small x, many of 
us believe that the dynamics that leads to unitarization can be understood 
in weak coupling [23 . It is interesting to speculate whether the physics of 
unitarization via weak coupling may also provide insight into the success- 
ful phenomenology of high energy scattering in the language of Pomerons 
and Reggeons. Historically, this physics was believed to be entirely non- 
perturbative in nature. 



ZEUS+Hl 




lo"** 10'^ 10'^ 10 ■' 



Figure 2: HI and ZEUS QCD fits to F2 data for fixed as a function of x. 

DIS experiments at HERA ushered in a new era in high energy QCD 
when they showed conclusively that structure functions rise very rapidly at 
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small X. It came as a surprise (even though it shouldn't have), since predic- 
tions for this behavior existed right from the very early days of QCD j24j . 
It has since been shown that DGLAP based QCD fits are very successful 
in fitting the HERA data |25j. These are shown in Fig. ^ In DGLAP evo- 
lution, the cross-section grows as one increases Q^. However, even though 
the multiplicity apparently grows with increasing Q^, the phase space den- 
sity actually decreases. This is because even though one sees many more 
partons, they are smaller in size. A cartoon demonstrating this is shown in 
Fig. ini We would like to ask what would happen if one were to fix and 
decrease x. In DGLAP, in part, this depends on the initial x distribution, 
which has to be a fairly steep power law. A natural mechanism is given by 
BFKL evolution, which gives rise to the rapid growth in x similar to those 
assumed in DGLAP initial conditions. The phase space density in this case 
grows very rapidly. Thus the great value of the BFKL equation is that it 
provides a mechanism to rapidly generate large numbers of partons, which 
in turn give rise to large cross-sections ^. 



A cartoon demonstrating this growth in the phase space density is shown 
in Fig. 01 When the phase space density becomes large, of order 1/as, par- 
tons fill up the entire phase space of the hadron saturating the phase space 
density and leading to the phenomenon known by this name [O] • Saturation 
is driven by higher twist effects corresponding to the screening and recombi- 
nation of partons in the high density environment. These contributions are 
non-linear in the parton density and their net effect is to slow down the rise 
in the parton density 0. These effects are not contained in BFKL evolution, 
which describes only linear evolution in the (unintegrated) parton densities. 

^Our perspective of BFKL here is in the S-channel dipole scattering picture of 
Mueller ^ 





Figure 3: Phase space density from DGLAP evolution. 
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The competition between the perturbative Bremsstrahlung and the many 
body, non-hnear screening and recombination effects is characterized by a 
scale Qs{x) in the x — plane. When Q < Qs{x) one expects the latter 
to dominate, while for Q » Qg, the linear evolution characteristic of QCD 
Bremsstrahlung is apparent. Saturation therefore provides a natural mech- 
anism to unitarize cross-sections. (The converse is not necessarily true: a 
unitarized cross-section is not necessarily a saturated cross-section.) 



The Color Glass Condensate is an effective field theory which describes 
the phenomenon of saturation in QCD |H| EI- The large parton density 
provides a large scale, the saturation scale oc UgNcdN / dy / tiE? , which 
ensures that the physics can be described in weak coupling. The high oc- 
cupation number, n ~ l/as, ensures that the small x modes can be treated 
as classical fields. These are coupled to static sources at larger values of x. 
The evolution of the source density with x is described by renormalization 
group equations. The strong fields in the CGC description ensure saturation 
already at the classical level of the theory. This is preserved by quantum 
evolution, which introduces additional non-trivial features. We will discuss 
the CGC further in Lecture IL 

At what values of x does saturation set in? This is a difficult question to 
answer from first principles ^ since one does not have quantitative control 
over where the higher twist effects that contribute to saturation become im- 
portant. The data presented in Fig. |2l don't appear to show any saturation 
effects. However, there are indications from the behavior of the gluon dis- 
tribution (at small x and moderately small Q^) and more importantly, from 
the longitudinal structure function Fl extracted from the data that higher 

^It is no easier to answer this question than to state at what values of perturbative 
QCD is applicable. The answer depends on what process one is calculating. 
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Figure 4: Phase space density from BFKL evolution. 
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twist effects may be playing a role already at HERA. It has been argued 
by Bartels, Golec-Biernat and Peters [^Tj that higher twist contributions 
that may be significant in Fi and Ft separately, cancel in F2 {=Fl + Ft)- 
Unfortunately, the Fl data from HERA has large systematic errors at small 

X. 

It is well known that the virtual photon-proton cross-section at small x 
can be written as |13 EHl EBl 



^T,L = J ^ y dz\^pT,L{r±,z,Q )| aqqN{r±,x), (3) 

where IVt.lP is the probability for a longitudinally or transversely polarized 
virtual photon to split into a quark with momentum fraction z and an anti- 
quark with momentum fraction 1 — z of the longitudinal momentum of the 
virtual photon. This probability is known exactly and it is convoluted with 
the cross-section for the qq-pa.ir to scatter off the proton. A simple model 
of small x DIS was introduced by Golec-Biernat and Wusthoff ^{. In 
this model, the qqN cross-section (often called the dipole cross-section) was 
parametrized to be 



where 



1 — exp ' 



(4) 



q's{^) = qU^) . (5) 



A 

X J 

Golec-Biernat and Wusthoff found a fit to all the HERA data for x < 10^^ 
and < 20 GeV^ with the parameters Qo = 1 GeV, ao = 23 mbarns, 
xq = 3 ■ 10~^ and A = 0.3. Interestingly, this model also provided a good fit 
to the diffractive data .30; , in particular the surprising energy dependence of 
the diffractive cross-section. The form of the hadronic cross-section adopted 
by Golec-Biernat and Wusthoff arises naturally in the CGC picture [SJ. 
The latter, as we shall discuss further in Lecture II, also has the right 
large k± behavior, the absence of which is one of the shortcomings of the 
Golec-Biernat- Wusthoff model. An improved version of the Golec-Biernat- 
Wusthoff model was developed by Bartels, Golec-Biernat and Kowalski |32j . 
who matched the saturation model to DGLAP evolution at higher values of 
Q^, while maintaining its key features. A further improvement of the model 
was the inclusion of the impact parameter dependence of the saturation 
scale by Kowalski and Teaney [33j . An alternative approach was followed 
by Mueller, Munier and Stasto [HI], who inferred the S-matrix (and therefore 
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the saturation scale in the Golec-Biernat-Wusthoff parametrization) from 
the t-dependence of exclusive /?-meson production. The value of the satu- 
ration scale extracted at small x in all the approaches lies in the ball park 
of ~ 1-1.5 GeV^. A similar analysis of J/tp production was performed 
by Guzey et al. [33 • For other saturation approaches, see Ref. [361- For 
a critical discussion of different models, see Ref. [57] . These authors con- 
clude that, for the most central impact parameters, the parton densities are 
large enough at HERA to generate semi-hard scales which significantly im- 
pact the dynamics. An extensive recent theoretical review of vector meson 
production of HERA can be found in Ref. j38| . 

The Golec-Biernat model also inspired Golec-Biernat, Kwiecinski and 
Stasto 0nj to plot the HERA data for the virtual photon-proton cross- 
section in terms of the dimensionless scaling variable r = Q^/Q^. Their 
result, as shown in Fig. [51 demonstrated that the data for x < 0.01 and 
0.045 < Q2 < GeV^ scaled very nicely function of r. Why these 
data scale up to these high at small x will be discussed in Lecture HI. 

f .0-' 




Figure 5: Geometrical scaling of the virtual photon-proton cross-section 
with T = Q^Ql 

A key ingredient in the Golec-Biernat-Wusthoff model is the dipole cross- 
section aqqp in Eq. |31 for which they introduced the simple model in Eq. ^ 
We will see in the following that the dipole cross-section is ubiquitous and is 
a common ingredient in both DIS and hadronic scattering at high energies. 

3 Lecture II 

In this lecture, we shall discuss an effective field theory for high energy 
scattering-the Color Glass Condensate. A key ingredient in our discussion 
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will be the n-gluon component of the hadron wave function (where n is 
large!). In order to construct this wave function, we need to quantize the 
theory on the light cone. The initial value quantum problem is formulated on 
an equal light cone time surface x+ = (t + z)/y/2 = 0. The great advantage 
of this approach is that the vacuum simplifies greatly ^ . For instance, the 
light cone quantized boost operator commutes with the light cone QCD 
Hamiltonian. This property is not satisfied by the usual equal time boost 
operator, which leads to particle creation under boosts. Thus if we wish to 
construct boost invariant wave functions, we must need do so on the light 
front. 

Light cone quantized quantum field theories have remarkable proper- 
ties It was first noticed by Weinberg j42j that the Feynman rules of 
scalar field theories simplified greatly in the limit where the light cone mo- 
mentum oo. Subsequently, Susskind showed that there exists an 
exact isomorphism between the Poincare group of a quantum field theory 
on the light cone and the Galilean subgroup of two dimensional quantum 
mechanics. Thus for instance, as in quantum mechanics, the QCD light cone 
Hamiltonian (-PqcD' generator of translations in x"*"), for instance, can 
be written as Pqqj) = -Pqcd + ^QCD- One can construct a complete set 
of Fock states that are eigenstates of the free Hamiltonian. Now since the 
light cone vacuum is an eigenstate of both -Pqqd and Pqq^, one can write 
the physical eigenstate in terms of a complete basis of "parton" eigenstates 
of bare quanta. This simple observation (coupled with the physics of time 
dilation) is the basis of the parton model in quantum field theory. An 
elegant paper by Bjorken, Kogut and Soper |3S] which demonstrates these 
features in Quantum Electrodynamics on the light cone is recommended 
reading for anyone interested in understanding high energy scattering in a 
light cone quantized gauge theory. 

In this light cone framework, the wavefunction of a high energy hadron, 
can therefore be expressed as 

\h >= \qqq > +\qqqg > H h \qqqg ■ ■ ■ qqggg > , (6) 

Each wee parton in a configuration containing a large number of partons 
carries a small fraction x = k~^/P~^ of the total momentum P+ of the 
hadron. These small x configurations (shown in Fig. ^ forming a fuzzy, 
strongly correlated parton cloud extending beyond the Lorentz contracted 

®As the reader might suspect, there is no free lunch and the complications of the 
usual vacuum are transferred elsewhere. However for the purposes of perturbation theory, 
perhaps to our eventual peril, these may be neglected. 
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Figure 6: Wee partons in the infinite momentum frame. 

valence partons, can only be probed at high energies, and are the relevant 
configurations for multi-particle production. We cannot at present write 
down the wave function for these configurations in a hadron ^. One can 
instead write down a path integral for the ground state of a hadron, which 
picks out the contributions from multi-parton states. 




Figure 7: A large nucleus boosted to high energies. A colored probe with 
transverse resolution x_l ~ << j^^^ sees a large number of color charges. 

A classical effective theory for the small x modes of a large nucleus was 
constructed by Larry McLerran and myself [Hj. This theory is a coarse 
grained field theory constructed in the infinite momentum frame (IMF) 

— > oo and in the light cone gauge A'^ = 0. In the IMF frame, only one 
component J"*" of the valence parton current is important-the others are sup- 

^However, this is precisely what Al Mueller has done for the specialized case of the 
n-gluon component of a large onium pair in the large Nc limit |2(i| . 
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pressed by 1/P~^ . Also, as shown in the cartoon in Fig.[7| due to Lorentz con- 
traction, a wee parton with a transverse resolution x± ~ l/k± « 1/Aqcd 
will see a large number of color charges from nucleons, all squished together 
in the transverse plane. 




Figure 8: The valence, glue and sea distributions in a proton. Note that 
the glue distribution is scaled down by a factor of 20! 

Another important ingredient in our construction is the observation that 
the valence and glue (+ accompanying sea) modes are cleanly separated in 
X. This is seen in Fig. |5J A remarkable feature of this plot is the extent 
to which the gluon distribution dominates the valence distributions as one 
goes to small x. Indeed, even at x ~ 0.1 there are twice as many gluons as 
up and down valence quarks combined. Therefore at small x we are con- 
structing a theory where the dynamical degrees of freedom are gluons. The 
valence partons act as sources and the following simple kinematics suffices 
to convince one that they are static sources: 

1 _2k+ _ 2xP+ 

n, J_ I 

2P+ _ 

,2 ^ ^wee ''"valence- {') 

Not much changes with the valence parton distributions over the time scales 
of interest for the dynamics of the wee partons. Thus one has something 
akin to the Born-Oppenheimer approximation-an important ingredient in 
an effective theory. However, one cannot integrate out the valence sources 



''"wee 



'^valence 
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completely out of the theory. This is because they carry color charge (lots 
of it as we shall see) and therefore couple to the wee partons ®. 



color charges 
at large x 



small X gluon 



•V • • • 

• >• • • • • 

• • • • 
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Figure 9: A long wavelength wee parton sees a lot of static charges (cx ^^/^) 
when it resolves small transverse distances. 

Since the wee parton has a large wavelength in the longitudinal direction, 
it can resolve a lot of color charges provided its transverse wavelength is not 
too large. The inequality 

^wee — Avalence — , \0) 

suggests that wee partons with x « A~^/^ can resolve partons all along 
the longitudinal extent of the nucleus. This is shown in Fig. El Here nip 
is the nucleon mass and 7 ~ /rriN is the Lorentz factor in the infinite 
momentum frame. However, if the wee parton had a wavelength k±_ < 
^QCD ~ 1 fill) it would see no color charge at all since color is confined (in 
nucleons!) on this scale. It is only if the wee parton has a short wavelength 
in the transverse direction k± >> Aqcd that it will see color charges from 
different nucleons along the longitudinal direction. These charges will be 
random since they are confined to different nucleons and do not know about 
each other. 

How many of these random sources the wee partons actually couple 
to depends on the typical transverse momentum of the wee parton ^. A 
wee parton with momentum p± resolves an area in the transverse plane 



We note here that the phenomenon of "Umiting fragmentation" observed in rapidity 
distributions over a wide swath of energies and reactions is basically a statement about 
the recoil-less universal nature of the valence ("fragmentation") distributions. 

^The wee parton is soft only in its longitudinal momentum-its transverse momentum 
may be large. 
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(Axj_)^ ~ The number of valence partons it interacts simultaneously 

with is then 

k = k^Ax^ f = ^'j^a'''' i^x^)"^ , (9) 

which indeed is proportional to A^^^ since A'^aience = 3 • j4 in QCD. 

For a large nucleus with k » 1, one can show that the most likely 
representation is a higher dimensional classical representation of order ^/k. 
The net charge that the wee partons couple to is a classical color charge and 
is represented by a classical color charge density (per unit transverse area) 
we denote as p. This argument can be made rigorous for SU{Nc) for large 
numbers of random sources 05] • Since the charges are random, we have 

< >= ; < p''{x^)p\y^) >= p\6-'6^'\x^ - y^) , (10) 

where 

""'hp- <"> 

is the color charge squared per unit area. For a very large nucleus, j4 >> 1, 
oc A^/^ » Aq(2;j-, is a large scale. Since it is the only scale in the effective 
theory, we expect that as{n\) « 1. Thus in the large A limit, one can 
construct the small x limit of QCD as a weakly coupled effective field theory. 

After these preliminaries, we can now write down the generating func- 
tional for the small x effective action, 

m = / [dp] w^, [p] J X /\ , , (12) 

J [ P [dA]6{A+)e''^l^'P^ J 

where A"*" is the longitudinal momentum scale separating the sources from 
the fields. Wa+[p] is a gauge invariant functional describing the distribution 
of sources at the scale A"*". The small x effective action can be written in 
terms of the sources and the fields as 

S[A,p] = ^J + ^ J d^x±dx-5{x-)Tv {pU-^,oo[A-]) , 

(13) 

where 

f/-oo,oo = Vexp (ig [ dx+A-'''TA (14) 



is a path ordered exponential along the light cone time direction. The second 
term in the effective action can also be written |17j as Tr (/9ln([/_oo,oo))- To 
the order studied thus far, the two forms of this term, which represents 
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the coupling between small x fields and large x sources, are equivalent. It 
remains an open question whether the latter form leads to different results 
at higher orders. 

As discussed previously, for a large nucleus, the sources are random light 
cone sources. For SU{Nc), this represents a random walk in the space of 
Nc — 1 Casimirs. It can be shown explicitly that the quadratic Casimir 
dominates the random walk (with the higher Casimirs giving contributions 
that are suppressed by powers ofl/Vk. The weight functional in the path 
integral, as first conjectured in Refs. [Hllini, can therefore be represented as 
a Gaussian functional weight (for large and one obtains, 

^[p]=exp(-y"dV^^^|^) , (15) 

where fi\ was defined in Eq. ^2 In general, W[p\ is not a Gaussian as we 
shall discuss shortly. 

Before we go there, let us first discuss the effective action in Eq . 111^1 The 
classical equations of motion for a fixed configuration of /?'s is given by the 
saddle point of the effective action in Eq. 1131 These are just the Yang- Mills 
equations, 

D^F^-'- = 6''+ 5{x-) p''{x±) . (16) 

The solution of the Yang-Mills equations are the non-Abelian analog of 
the Weizacker- Williams fields in classical electrodynamics 0. As one may 
recall, when one boosts a classical charge in electrodynamics to the IMF, 
the fields of the charge look like that of a sheet of plane polarized radiation. 
The fields are singular on the sheet and pure gauges outside. This is shown 
in Fig. Uni 

The solution is given by A~ = and 

All = J^'P exp (^ig:^p{x±,x')^ Vj exp (^ig:^p{x_L,x')^ . (17) 

The path ordering here is in x~ . Note also the x~ dependence in p. A 
careful solution |481 149j of the Yang-Mills equations requires smearing in 
x~ . A final note on the solution is that the charge density p that appears in 
the solution is not the color charge density in light cone gauge but instead 
the color charge density in covariant/Lorentz gauge. In the latter, one has 
the solution A'^ = ^p6{x-), A'~ = A! x. = 0. 

The explicit solution of the gauge field in terms of p is not trivial. How- 
ever, since we are interested in number distributions, one can simply replace 
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a_i (pure gauge) 



boost ( y= 3 ) 



Figure 10: Weizsacker- Williams fields of a charge boosted to the infinite 
momentum frame are pure gauges on either side of the charge. 

the measure [dp] [dp] in the path integral. The Jacobian in the trans- 
formation is quite simple '48 and does not contribute to the result. Thus 
one can compute distributions in light cone gauge quite straightforwardly 
by expressing them in terms of charges in covariant gauge. To be specific, 
one averages the solution in Eq. El with the weight functional 

< AA >p= I [dp]A,i,[p]A,i,[p]W^+[p] . (18) 




Figure 11: The occupation number of classical non-Abelian Weizsacker- 
Williams fields. 
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For the Gaussian weight in Eq. 1151 which is appropriate for large nu- 
clei and not too small x, one can analytically compute the solution. The 
result is shown in Fig. 1111 The occupation number is defined to he (j) = 
{2TTfdN/7rR'^/(Pk^dy/2{N^ - 1) with Ql « a^iVc/i^ ln(Q2/A2^^p). Hence 

A^/^ln{A) ~ A^^^ for A >> 1. At large transverse momenta {k± >> 
Qs), the distributions have the characteristic oc fj,\/k'j_ form of Weizsacker- 
Williams gluons in electrodynamics. At smaller transverse momenta {k± « 
Qs), the distribution has the form cp ~ ln(Qs/k±)/as. Thus in light cone 
gauge, the strongly non-linear behavior of the fields is responsible for the 
softening of the infrared behaviour of the classical fields. This non-linearity 
is responsible for the phenomenon of saturation. A full description however 
requires a discussion beyond the classical level discussed thus far. 

We are now in a position to understand the term Color Glass Condensate 
for the state we are describing, namely, the ground state properties of a 
hadron/nucleus at very high energies. Color is obvious since the state is 
comprised of a large number of gluons. It is a condensate because, as we 
have just seen explicitly, the gluons have occupation numbers (p ^ 1/as and 
have momenta peaked at k± ^ Qg. Finally, it is a glass because the gluons 
are coupled to random sources with time scales of evolution much longer 
than those of natural time scales associated with the scattering. 

We have discussed thus far a classical effective field theory for large nu- 
clei and Gaussian sources. Though many features of the theory persist, the 
theory is inadequate to describe small x evolution at very small x. The main 
culprit is the Gaussian assumption for W^[/9]. One does not obtain Gaussian 
correlations when quantum corrections are included. The first suggestion 
this was the case came from computing small quantum fiuctuations about 
the classical saddle point solution. It was noticed that the corrections, pro- 
portional to as ln(l/x), became arbitrarily large at small x [50]. It therefore 
was important to sum these up a la BFKL jlH]- The method that was in- 
vented to do this, now called by the acronym JIMWLK [5] after the names 
of the principal contributors, was a Wilsonian renormalization group ap- 
proach. Small fluctuation corrections to the effective action in Eq. ^Jat one 
step in X (A^) are incorporated into a modified source density at the next 
step: Wyv+[p] — > W^A'+I/^']' where p' = p + 5p, is the new classical source 
density which incorporates the small fiuctuations in the fields at the scale 
A"*" into the source density at the new lower momentum scale A'^. (One 
requires that 05 ln(A^/A''^) << 1 to ensure that the small fiuctuations are 
under control at each step.). These are shown in Fig |12l 

The JIMWLK evolution equation for the RG evolution of the weight 
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Figure 12: The weight functional modified under evolution. 



functional in Eq. El can be written as 



where X°'^{xL:yL)[p\ =< ^p"'{x±)6p^{y±) >p is a two point function in the 
background field of the hadron. From Eq. I19L one can construct a master 
equation for n-point correlators. With a change of variables, a° where 

V^a = p, one obtains for the expectation value < 0[a] >y of an operator 
O the relation, 

d < 0[a] >Y 1 /■ ^ abr ^ 1 rnn\ 



This equation is a generalized Fokker-Planck equation in functional space, 
where Y is "time" and x is the diffusion coefficient |511 I52j . x has the 
functional form. 



where 



f (x-z)-(y--z) ) ( t )(i_yty 

(21) 

V^{x^) =Ve^p(^ig J dx' a"" {x~ , x ±)T''^ , (22) 

and V denotes path ordering in x^ . 

Consider for instance the two point function < a{x±)a{y±) >y, we dis- 
cussed previously. In the weak field limit, ga « 1, the JIMWLK equation 
for this two point function reduces to the BFKL equation '21' . 

Interestingly, one can solve the JIMWLK equations for the two point 
functions in the other limit-the strong field limit-gfa ~ 1, using a mean 
field approximation, which we shall call the Random Phase Approximation 
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(RPA) ^2 . The phase space density (/> we discussed in the classical approx- 
imation, has the following approximate solutions at large rapidities Y: 



Here u; = 4 In 2 is a constant and as = asNc/TT. The phase space density 
in the first line (at very large k±) corresponds to the MV/DGLAP regime 
discussed previously. The regime described in the second line here is new 
and is the BFKL regime of quantum evolution. The last line corresponds to 
the dense Color Glass regime of k_i « Qg. Interestingly, the phase space 
density in this dense regime has the same structure as in the classical theory 
except that the color charges are screened at distances greater than 1/Qs [SSI 
l5l] . as opposed to Aqcd in the classical theory. A good approximation to 
the behavior of the phase space density in the three regions is the form |53j 



where varying the anomalous dimension 7 from 0.63 in the "BFKL" region 
to 1 in the MV/DGLAP region. We will discuss this point further shortly. 

The JIMWLK equations are master equations for n-point correlators 
under small x evolution. These do not have closed form expressions. For 
instance the evolution equation for 2-point correlators contains 4-point cor- 
relators, and so on. This hierarchy is analogous to the well known BBGKY 
hierarchy in statistical mechanics. Like the latter, one has to make as- 
sumptions about correlators at a certain order to close the hierarchy. The 
JIMWLK equations have not been solved analytically even though there has 
been an attempt to solve them numerically jl8j . 

The Balitsky-Kovchegov equation J2] is a simple non-linear evolution 
equation for the forward scattering amplitude. The derivation by Kovchegov 
was performed in the framework of Mueller's dipole picture |'26^. It is equiv- 
alent to the JIMWLK expression for the forward scattering amplitude in the 
limit of large Nc and a'^A^^^ >> 1. Recall the expression we had in Lecture 
I for the virtual photon-proton cross-section (Eq. The hadronic (dipole) 




(23) 




(24) 
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cross-section there can be written as |31| 

cTgqp{x,r^,b) = 2 j dH{l-ReS{x,r^,h)) , (25) 

where the S'-matrix can be written in terms of the path ordered exponential 
in Eq.lHlas 

ReS{x,r^,b) = ^ < Tr {V\x^)V{y^)) >Y= I - M{x,r^,b) , (26) 

where M is the imaginary part of the forward scattering amphtude, rj_ = 
x±—y±, b = {x±+y±)/2 and the path ordered exponential V was introduced 
in Eq. 1221 Balitsky and JIMWLK have shown that the S-matrix in Eq. I^Hl 
satisfies the equation 

QY - 27r y ^ (x - z)2(z - y)2 < ^'^^"^ ^-Mn^. Vy) -Lr^K, Vy) > . 

(27) 

Here = agNc/ir. 

In the Nc ^ 00 and a'gA^^^ 00 limits, the average over the product 
of traces here can be written as the product of the averages of the traces 

< TriV^V^mV^Vy) >=< Tt{V^V,) >< TriV^Vy) > . (28) 

With this factorization, one now obtains the Kovchegov "mean field" equa- 
tion for the imaginary part of the forward scattering amplitude 



dMix±,y±) _ f ix±-y±) 

— a 



2 

My{x±,Z±) +AAy(y_L,^±) - AAy(x_L,y±) 



dY J^^ (x_L - z±)'^{z± - yi_Y 

- AAy(xx,z^)-AAy(zx,y±)] . (29) 

In the limit where A/" << 1, the non-linear term in Eq. [551 can be ignored 
and the equation reduces to the BFKL equation. In this limit, the amplitude 
has the solution, 

A/V(.-l) » (riQji)'/^e"°'>-exp(- '"'^y:^^°' ) 



^°One would imagine that in QCD that one would need only the large limit for this 
factorization. This can be seen for instance from the one loop effective action in QCD |39) . 
However, in our case, we are dealing with a very particular background field which does 
not factorize for large Nc alone. It requires the large A limit as well for factorization to 
hold. I thank F. Gelis for an enlightening discussion on this point. 
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where w = 41n2 2.77, /? = 28C(3) ^ 33.67 and p = ln{r^Qof. If we define 
the saturation condition as 

the vanishing of the exponent for this value of r± gives 

A more careful solution of the BK-equation close to the saturation boundary 
gives c = 4.88. 

We can replace Qo in Eq. EEl by Qs from Eq. 021 One thus obtains 

^fy{r^) - {rlQlV^ e-^^ exp ("i^^^^^^) , (33) 
with 7s « 0.63. In the strong field RPA limit [521 ^DIM]) one finds 

ATy = 1 - Kexp (^-j- ln2(riQ2)^ ^ (34) 

with c = 4.84 and k is an undetermined constant. A complete analytical 
solution of the BK equation is still lacking though, as we shall discuss, there 
have been interesting recent developments in this direction |64| I65| I66j . 

The Balitsky-Kovchegov equation has however been solved numerically |13l 
I14 | I15 | ITH]. These solutions have the following features: 

• The dipole amplitude is shown to unitarize and solution is shown to 
exhibit geometrical scaling. 

• The saturation scale has the behavior predicted in Eq. 1321 

• The infrared diffusion problem of the BFKL solution is cured by the 
non- linear term in the BK equation. 

• Ultraviolet diffusion is still present in the BK equation. This can be 
cured by including running coupling effects jl4l I17j . 

Numerical solutions of the Balitsky-Kovchegov equation are shown in Figs. 1131 
and 1141 In Fig. 1131 the unintegrated gluon distribution is shown to be sta- 
ble against infrared diffusion. The unintegrated gluon distribution from the 
BFKL equation is shown for comparison. In Fig. 1141 the solution is shown 
explicitly to exhibit geometrical scaling for k± < Qg. The solutions also 
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Figure 13: The unintegrated gluon distribution (j) from numerical solu- 
tions ^3 of Balitsky-Kovchegov equation plotted as a function of trans- 
verse momentum. The dotted lines lines are the BFKL results while the 
solid ones are from solutions of the BK-equation. 




k/Q,(Y) 



Figure 14: Geometrical scaling of the unintegrated gluon distribution cf) with 
Qs/k from numerical solutions |E] of the Balitsky-Kovchegov equation. 
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exhibit a soliton-like structure; as shown by Munier and Peschanski, this 
analogy is an exact one . Numerical solutions to the full JIMWLK equa- 
tions have been studied by Rummukainen and Weigert Much work 
remains in that direction. 

As we saw in Fig. |3 the HERA data appears to exhibit geometrical 
scaling up to rather large values of Q^, values significantly larger than Q^. 
Why is this so? A qualitative explanation was provided by lancu, Itakura 
and McLerran (See also Ref. j59j.^ We can write the solution of the 

BFKL equation in Eq. EOl as 

A/,.(r.)=»exp(.5.y-f-^), (35) 

where p = ln{l / r'j_QQ) . Now from our definition of the saturation scale 
(Eq. 021), we can write p = pc + asY, where pc = ln(l/r^(5g). One finds 
that the solution scales as long as /j^ << 2f3asY. Since we are interested 
in the region » Q^, we find that geometric scaling holds for XasY « 
In(QVQo) << ^oisY + ^/2l3asY. This is clearly valid for V2jMFY » 1 
and is thus easily satisfied even for moderate Y^s. 

A quantitative fit to HERA data using a dipole parametrization match- 
ing the BFKL dipole with the mean field (RPA) form in Eq. |^ was per- 
formed by lancu, Itakura and Munier [HSl- They obtain good fits but the 
HERA inclusive data at present cannot distinguish between these fits and 
models that do not contain saturation [^. In general, a large number 
of models give good fits to the inclusive data. Its the non-inclusive data 
that provides a more sensitive test and preliminary evidence there is that 
the CGC models fare better than those without saturation effects jSZl- To 
summarize, as discussed in the previous lecture, CGC based models do rea- 
sonably well but with the cautionary caveat that much more work remains 
for a consistent phenomenology. 

How does Qs behave as a function of rapidity (or energy)? For fixed 
coupling, as we discussed previously, we obtain Eq. 02 If we allow the 
coupling to run, we get very simply, 

Qlvnnnmgas = ^QCD exp (^^/2boc{Y + Yq)^ , (36) 

where bo is the coefficient of the logarithm in the one loop QCD /5-function. 
The former expression can be understood as the low energy (rapidity) limit 

^^The numerical solutions of course provide a quantitative answer, but don't necessarily 
provide insight into the phenomenon. 
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of the latter. Now, the inclusion of running coupling effects in BFKL is 
highly non-trivial. The proper treatment of next-to-leading order BFKL re- 
quires renormalization group improvement of the NLO BFKL kernel, where 
all collinear singularities are summed to all orders j^. It was suggested 
by Mueller and Triantafyllopolous that the saturation scale could be ex- 
tracted from solving the BFKL equation in the presence of an absorptive 
boundary [201 • Triantafyllopolous showed that this treatment could be 
extended to the resummed NLO BFKL case, and he was able to extract the 
energy dependence of Qs- His result is shown in Fig. E| The value of A is 
rather large for the fixed coupling case. Furthermore, pre-asymptotic effects 
are clearly seen to be important. They are less so for running coupling. The 
"naive" running coupling result is seen to be a fairly good approximation to 
the full resummed result. The latter is seen to be to be quite close to the 
value extracted from the HERA experiments. 



The Logarithmic Derivative of Qs^ 
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Figure 15: The energy dependence of the saturation scale expressed in 
terms of A = d\nQl/dY . 

The A-dependence of the saturation scale Qs has been computed recently 
by Mueller in this picture . The result is shown schematically in Fig. 1161 

This result has the form « ^qcd e^^^^^^ {A^/^) ^ where ^ here is an 
undetermined constant. At small rapidities, the second term dominates and 
one has the result we discussed previously in lecture H; namely, cx A^/^ . 
For large Y , the converse is true, and the A dependence of the saturation 
scale is gradually lost. Thus as illustrated in Fig. El for large rapidities, one 
may find the remarkable result that the saturation scale, for fixed impact 
parameter, is independent of the target. 

Finally, we now turn to very recent developments. As we mentioned 
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universal fixed point? 




sqrt(Y) 

Figure 16: The A-dependence of the saturation scale as a function of the 
square root of the rapidity. The initial relative height is of order A^^^. 

previously, Munier and Peschanski made the remarkable observation that 
the BaHtsky-Kovchegov equation, in a diffusion approximation (see footnote 
3), can be shown to he (with appropriate field re-definitions) in the same 
universality class as the Fischer-Kolmogorov-Petrovsky-Piscunov (FKPP) 
equation inZl; which describes the behavior of travelling wave fronts. So- 
lutions to the latter are known, and Munier and Peschanski were able to 
deduce from these the asymptotic form of the amplitude in both the fixed 
and running coupling B-K equation as well as the form of the saturation 
scale. They were also able to compute the first two universal pre-asymptotic 
contributions (in rapidity) to the saturation scale. The first pre-asymptotic 
correction was also computed by Mueller and Triantafyllopoulos [3^1 IHUj . 
whose approach, as we discussed previously was apparently completely dif- 
ferent. 

Thus far, we have not discussed the impact parameter dependence of 
the amplitudes. As we shall see, this ties into the previous discussion in 
an interesting way. It was observed by Mueller and Shoshi |2H] that, for a 
fixed impact parameter, the amplitude AA(?)) could violate unitarity (in dense 
"hot spots"), even if the impact parameter averaged amplitudes remained 
well below unity. The possibility that this could happen in dipole-dipole 
scattering was suggested by Mueller and Salam 1^2] and observed in numer- 
ical simulations of the same by Salam [TDj. Mueller and Shoshi suggested 
that these hot spot violations of unitarity could be corrected for by solving 
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the BFKL equation with a second absorptive boundary at AA = a| in 
addition to the first one at ~ 1. Their computation of the saturation 
scale showed large corrections which vanished only at asymptotically small 
values oi as- 

It was noticed subsequently by lancu, Mueller and Munier 6'5_ that 
these results could be interpreted in terms of a generalization of the FKPP 
equation, the stochastic FKPP (sFKPP) equation, well known to experts in 
statistical mechanics studying the properties of traveling wave fronts j71j . 
It was realized in these works that the original FKPP equation, in the re- 
gion where the wavefront is small, did not properly treat the contribution of 
diffusion to the evolution of wavefronts. It is this diffusive property which 
provides the seed for the evolution of wavefronts in regimes where the oc- 
cupation number is extremely small [Z2l- This diffusion is modeled by a 
stochastic Langevin term in the KFPP equation. The picture described by 
this sFKPP equation is now that one has an ensemble of wavefronts, each 
one stable, and saturation scales that are Gaussian distributed among the 
fronts. A very recent comprehensive analysis of the implications of these 
ideas for BK and JIMWLK has been provided by lancu and Triantafyl- 
lopoulos [nni- These are severe. While the physics in the saturation regime 
at high energies is likely well described by BK and JIMWLK, the compe- 
tition between the dispersion of wavefronts (given by the variance of the 
Gaussian) and the scale governing the exponential decay of fronts (the BK 
anomalous dimension) suggests that there may be no BFKL regime at any 
energy. lancu and Triantafyllopoulos have suggested a Langevin generaliza- 
tion of Balitsky's hierarchy which includes these effects. The ramifications 
of these ideas (particularly how they arise in the Color Glass Condensate) 
remain to be fleshed out further. However it should be clear to the reader 
that the stochastic nature of high energy processes and the exact analogies 
to well plumbed models in statistical physics is an exciting direction for 
research in high energy QCD. 

4 Lecture III 

We now come to the final topic of these lectures, high energy hadronic scat- 
tering. Consider for instance the scattering of two nuclei at high energies-a 
problem very relevant at high energies. In the CGC picture, an observable 

^^The second boundary also had the virtue of correcting for the frame dependence of 
the BK-equation at high energies. 
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< O >Y can be computed as 

<0>Y= I [dpi] [dp2] W^M WM 0{pi,p2) , (37) 

where Y = lii(l/xi7') and xp = xi — X2- All operators at small x can be 
computed in the background classical field of the nucleus at small x. In par- 
ticular, the Yang-Mills equations are solved to determine the dependence 
of the background field in terms of the sources All quantum infor- 
mation, to leading logarithms in x, is contained in the source functionals 

To be specific, inclusive gluon production in the CGC is computed by 
solving the Yang-Mills equations [D^,F'^'^]'^ = J'^'", where 

r = ppi6{x-)5^+ + Pp2Kx^W~ ■ (38) 

with initial conditions given by the Yang-Mills fields of the two nuclei before 
the collision. These are obtained self-consistently by matching the solutions 
of the Yang-Mills equations on the light cone. The initial conditions are 
determined by requiring that singular terms in the matching vanish. Since 
we have argued in the previous lectures that we can compute the Yang- Mills 
fields in the nuclei before the collision, the classical problem is in principle 
completely solvable. Similarly, quark pair production at small x can be 
computed by solving the Dirac equation (again by matching on the light 
cone) in the background field of the two nuclei with the latter obtained from 
solutions of the Yang- Mills equations. 

This approach therefore solves the problem of the behavior of wee partons 
first outlined by Bjorken nearly 30 years ago /75/. As we shall discuss in 
the following, this enables one to calculate from first principles the initial 
energy density and the formation time of partons in a high energy heavy ion 
collision (first estimated by Bjorken f76] /j. 

Hadronic scattering in the CGC is studied in practice through a system- 
atic power counting in the density of sources in powers of pi, 2/^1.1 2- This 
power counting in fact is more relevant at high energies than whether the 
incoming projectile is a hadron or a nucleus. In addition, one can begin 
to study the applicability of both collinear and k± factorization at small x 
in this approach. We shall discuss here how to interpret gluon and quark 
production at high energies in hadronic collisions for the cases when both 

^^For more detailed discussions of computing amplitudes in strong background fields, 
see Refs. [ZS1II1EI1IH3- 
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hadrons provide either dilute or dense sources for the scattering as well as 
the asymmetrical case of one source being dilute and the other being dense. 

4.1 Gluon and quark production in the dilute/pp regime: 

iPpi/klpp2/kl « 1) 

The power counting here is applicable either to a proton at small x, or 
to a nucleus (whose parton density at high energies is enhanced by ^4^/^) 
at large transverse momenta. The relevant quantity here is Qs, which, as 
one may recall, is enhanced both for large A and small x. So as long as 
k± » Qs >> Aqcd, one can consider the proton or nucleus as being 
dilute. 

To lowest order in ppi/k\ and Pp2/k\, one can compute inclusive gluon 
production analytically. This was first done in the = gauge j77j 
and subsequently in the Lorentz gauge d^A^ = |2Bj- At large transverse 
momenta, Qs « the scattering can be expressed in a A;_|_-factorized 
form shown in Fig. I17L The inclusive cross-section can be expressed as the 
product of two unintegrated {k± dependent) distributions times the matrix 
element for the scattering. This type of factorization (as distinguished from 
the collinear factorization of perturbative QCD) is called ^^-factorization. 
In this case, all the small x evolution can be factorized into the unintegrated 
gluon distributions. 

The comparison of this result to the collinear pQCD gg gg process 
and the k± factorized gg ^ g was performed in Ref. [73j. At this order, the 
result is equivalent to the perturbative QCD result first derived by Gunion 
and Bertsch [HHl- This result for gluon production is substantially modified, 
as we shall discuss shortly, by high parton density effects in the nuclei. 




Figure 17: factorized monojet. 
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One can show that k± factorization is a good assumption at large mo- 
menta for quark pair-production as well. This was worked out in the CGC 
approach by Frangois Gelis and myself [SJ. The result for the gauge field 
obtained in Ref. JTS, is put to use here. The result for inclusive quark pair 
production can be expressed in k± factorized form as 

dai 1 f (fki± (fk2± 



f d^ki±d/^ 

J (2vr)2 ^^^y S{k^^ + k2^-P±-1±) 



dypdyqdPp^<Pqi_ {2Tif Cl J (27r)2 (27r)2 

{\'n^ah{ki,k2]q,p) 



'fi{kn_)ip2{k2±) -^^—^ , (39) 



where and <j}2 are the unintegrated gluon distributions in the projectile 
and target respectively (with the gluon distribution defined as xG{x, Q"^) = 

d{k'j_) (l){x,k±)). The matrix element Tr ^\m~f^{ki,k2;q,p)\'^^ is iden- 
tical to the result derived in the ^^-factorization approach |82[ I83j . This 
result has been applied extensively to study heavy quark production at col- 

- - Tr (|m;^+(A;i,A;2;g,p)|^) 

lider energies In the limit \ki±\ , |A:2±| — > 0, — — , ^ is well 

defined-after integration over the azimuthal angles in Eg. 1391 one obtains the 
usual matrix element \-M-\'^g^gq, recovering the lowest order pQCD collinear 
factorization result. 



4.2 Gluon and quark production in the semi-dense/pA re- 
gion (pp/kl « IpA/kl ~ 1). 

The power counting here is best applicable to asymmetric systems such as 
proton- nucleus collisions, which naturally satisfies the power counting for a 
wide range of energies. Of course, as one goes to extremely high energies, 
it is conceivable that the parton density locally in the proton can become 
comparable to that in the nucleus. We shall discuss how that case works as 
weU 11 

In the semi-dense/pA case, one again solves the Yang-Mills equations 
[Df^, F^^'^] = J^, now with the hght cone sources J^'°- {=6^^ Pp{x±) + 

5'^^ 5(x+) p\{x_[_)), to determine the gluon field produced-to lowest order in 
the proton source density and to all orders in the nuclear source density. 
The computations are performed in Lorentz/covariant gauge dmuA^^ = 0. 

^*The reader might wonder about the frame dependence in this formulation. Our for- 
mahsm is vahd as long as the infinite momentum frame is applicable for both projectile 
and target, namely their momenta are much larger than the momenta of the constituents 
in the hard scattering. 
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The inclusive gluon production cross-section, in this framework, was first 
computed by Kovchegov and Mueller [S^ and shown to be k± factorizable 
in Ref. jBHl- In Ref. the gluon field produced in pA collisions was 

computed explicitly in Lorentz gauge. One obtains, 

with k2 = q — ki and U &z. V are path ordered Wilson lines containing all 
orders in the nuclear source density pA- The coefficient functions Cfj and 
Cv are simply related to the well known Lipatov effective vertex through 
the relation C7^ = + 

The path ordered exponentials U are color matrices arising from the 
rotation of the color charge density of the proton source due to multiple 
scattering off the nucleus. The path ordered exponentials V (differing from 
the C/'s by a 1/2 factor in the argument of the exponential) arise from the 
Green's function solutions of the equations of motion. Interestingly, they 
do not appear in the final result for gluon production. This is because for 
gluons produced on shell one finds remarkably that Cu ■ Cy = Cy = and 

= Cl = Akf^k'^^/q'^. Thus only bi-linears of the Wilson line U survive 
in the squared amplitude that gives us the gluon production cross-section. 
The result is /c^-factorizable analogous to Eq. except now one replaces 
02 with the unintegrated nuclear gluon distribution (pA ot< >p^. This 
distribution contains powers of the usual unintegrated gluon distribution 
to all orders-one recovers the usual unintegrated gluon distribution {(p2 in 
Eq. 139(1 at large transverse momentum. 

Our result in Lorentz gauge is exactly equivalent to that of Dumitru & 
McLerran in A'^ = gauge [89]. We now turn to a discussion of the Cronin 
effect which can be interpreted in this framework. 

The Cronin effect 

The Cronin effect was discovered in proton-nucleus collisions in the late 
70's |9U1 1911 I92j . The effect observed was a hardening of the transverse 
momentum spectrum in proton- nucleus collisions, relative to proton-proton 
collisions, that sets in at transverse momenta of order A:^ ~ 1 — 2 GeV, and 
disappears at much larger kj_''s. A corresponding depletion was seen at low 
transverse momenta, accompanied by a softening of the spectrum. At that 
time, the effect was interpreted as arising from the multiple scatterings of 
partons from the proton off partons from the nucleus [HSI- At high k±, the 
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higher twist effects, which, in the language of perturbative QCD, are respon- 
sible for multiple scattering [SH IHS] are suppressed by powers of k±. The 
relative enhancement of the cross-sections at moderate k±'s should thus die 
away ~ and indeed, the data seemed to suggest as much. Though a qualita- 
tive understanding of the previously observed Cronin effect was suggested by 
perturbative QCD, a quantitative agreement for all its features (such as, for 
instance, the flavor dependence) is still lacking. The Cronin effect has been 
observed at central rapidity at RHIC, and has been interpreted in terms of 
multiple scatterings, both in colhnear factorized pQCD [HHl 1^71 imn] 
and in saturation inspired models |1U11 11U21 11U51 11U61 ITB] . 
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Figure 18: Depletion of the Cronin peak from = to = 3 for minimum 
bias events. From Ref. jlU8j . 
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Figure 19: Centrality dependence of the Cronin ratio as a function of ra- 
pidity. From Ref. [TUH] . 

First data from RHIC on forward D-Au scattering at ^/s = 200 GeV/nucleon 
demonstrate how the Cronin effect is modified with energy or, equivalently. 
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with the rapidity. The x values in nuclei probed in these experiments, at 
A;_L ~ 2 GeV, range from 10~^ in the central rapidity region down to 10~^ at 
very forward rapidities A dramatic result obtained by the BRAHMS ex- 
periment (108) which has taken data up to pseudo-rapidities ri = 3.2 16 is the 
shrinking of the Cronin peak rapidly with rapidity and at higher rapidity, 
one sees in Fig. ^Jthat there is a significant suppression instead. Equally 
interesting is the centrality dependence of the effect. At central rapidities, 
one observes that the Cronin peak is enhanced in more central collisions, 
while, for forward rapidities, the trend is reversed: more central collisions 
at forward rapidities show a greater suppression than less central collisions! 
This is shown in Fig. 1191 

These results have a natural qualitative understanding in the CGC frame- 
work. Comparisons of saturation inspired models to the RHIC D-Au data 
lead to the following key conclusions: 

• Tree level partonic re-scattering responsible for the classical Cronin 
effect in high energy Deuteron-Gold collisions is computable in the 
CGC approach. When the weight for the color sources is Gaussian 
(which is valid at moderately small x for large nuclei), the formalism 
reduces to the familiar Glauber formalism of independent multiple 
scatterings. In the McLerran-Venugopalan model, a peak appears at 
k± ~ Qs,A 1 which is more pronounced for more central collisions [IDH 

[Tn^[Tie[M)| . 

• Quantum corrections (due to small-x evolution) to the tree level par- 
tonic re-scattering can be "naively" included by letting Qs — > Qs{x)- 
This model of quantum evolution predicts a larger Cronin effect at 
larger rapidities (because Qs grows as x decreases) sharply disagree- 
ing with the RHIC d-Au data at forward rapidities. On theoretical 
grounds it has been known for some time that quantum evolution 
quickly destroys ^'^ the Gaussian weight functional |5nj . 

• Proper quantum evolution in the CGC is described by the JIMWLK 

^^Guzey, Strikman and Vogelsang have recently argued that the average x of the 
BRAHMS data is significantly higher |1()7) . This is indeed true for the 2 — > 2 process. 
However, at small x, the 2 — > 1 process becomes important and samples smaller a;'s in the 
nucleus than in the 2 — > 2 process. 

^®The trends seen by BRAHMS are also corroborated by PHOBOS, PHENIX and STAR 
in different kinematic ranges |109l 111(11 

'^^This also explains why at sufficiently small x it would be more consistent to use 
something like the BFKL+RPA fit of Ref. [55] in interpreting DIS data than Golec- 
Biernat-Wusthoff type fits. 
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evolution equation. All simplified versions of JIMWLK (be they BK 
or RPA) lead to a "Cronin suppression" at large rapidities jll21 
I1U51 1161 1113j . Furthermore, the anomalous dimensions governing the 
hardness of the momentum spectrum discussed previously, of 

BFKL-type. These lead to a softer spectrum whose magnitude is more 
severely suppressed for more central collisions, exactly as predicted 
by the data. Another corroborative piece of evidence is the prelimi- 
nary observation of the broadening of azimuthal correlations between 
forward going and central hadrons in D-Au collisions by STAR |114j . 
This observation confirms a prediction by Kharzeev, Levin and McLer- 
ran |115j . 

From the qualitative features of the RHIC d-Au data the following pic- 
ture emerges The MV model works reasonably well at x ~ 10~^ - central 
rapidities at RHIC - where small x evolution is not significant. It is a good 
model of the initial conditions for quantum evolution as a function of ra- 
pidity. Correspondingly, the fact that the MV model fails badly at forward 
rapidities suggests that quantum evolution effects cannot be accounted for 
simply by a rescaling of the saturation momentum. For instance, calcula- 
tions with the RPA model reveal a qualitatively different behavior of the 
ratio RdA- In particular, a prediction of this mean- field solution is that the 
ratio RpA tends to A^"'~^^^^ at large transverse momentum, where 7 is the 
BK anomalous dimension. 

An important test of the CGC in D-Au collisions is provided by electro- 
magnetic probes |1U21 11U41 11161 1117j . If one observes the Cronin effect and 
its subsequent depletion with rapidity, one cannot but conclude that it is an 
initial state effect. Else, it may be predominantly due to final state effects, 
as has been proposed recently jll8j . 

Quark production in p/D-A collisions 

Quark production can now be computed with the gauge field in Eg.EOligS]. 
The field is decomposed into the sum of 'regular' terms and 'singular' terms; 
the latter contain 5{x~^). The regular terms are the cases where a) a gluon 
from the proton interacts with the nucleus and produces a qq-paii outside, 
b) the gluon produces the pair which then scatters off the nucleus. Naively, 
these would appear to be the only possibilities in the high energy limit where 
the nucleus is a Lorentz contracted pancake. However, in the Lorentz gauge, 
one has terms identified with the singular terms in the gauge field which cor- 

^^This picture will have to be substantiated by more detailed (read quantitative) com- 
putations and a thorough confrontation to the experimental data. 
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respond to the case where the quark pair is both produced and re-scatters 
in the nucleus! Indeed, the contribution of this term to the amphtude can- 
cels the contribution of the term proportional to the V's (see Eq. l4Up in the 
regular terms. 

Our result for quark pair production unlike gluon production, is not 
strictly k± factorizable. It can however still be written in k± factorized form 
as a product of the unintegrated gluon distribution in the proton times a sum 
of terms with three unintegrated distributions, (j)g,g, 4'qq,g and (pgqm- These 
are respectively proportional to 2-point, 3-point and 4-point correlators of 
the Wilson lines we discussed previously. For instance, the distribution (pqq^g 
can be interpreted as the probability of having a qq pair in the amplitude and 
a gluon in the complex conjugate amplitude. For large transverse momenta 
or large mass pairs, the 3-point and 4-point distributions collapse to the 
unintegrated gluon distribution, and we recover the previously discussed 
A;_L-factorized result for pair production fEq. l39|l in the dilute/pp-limit. Due 
to the length of the expressions, we will not present them here but refer the 
reader to Ref. 88 . 

Single quark distributions are straightforwardly obtained. Here the 4- 
point correlator (pqq^qq collapses upon integration over the momentum of the 
quark or anti-quark to the 2-point correlator (pq^q corresponding to a quark 
(or anti-quark) in the amplitude and complex conjugate amplitude. 

For Gaussian sources, as in the MV-model, these 2-, 3- and 4-point func- 
tions can be computed exactly as discussed in Ref. [SB] . Single quark distri- 
butions in the MV-model were recently computed by Tuchin |119) . Explicit 
computations of the size of /c^-violating terms as a function of quark mass 
and Qs are underway and will be reported shortly |124j . 

The results for gluon and quark production in Proton/ Deuteron- Gold 
collisions, coupled with the previous results for inclusive and diffractive \10Sl 
\125\. V^'A \12f)^ distributions in DIS suggest an important new paradigm. As 
one goes to smaller values of x in DIS and hadron colliders, previously 
interesting observables such as gluon distributions are no longer the right 
observables to capture the relevant physics. Instead they should be replaced 
by these dipole and multipole correlators of Wilson lines that seem ubiquitous 
in all high energy processes and are similarly gauge invariant and process 
independent. The renormalization group running of these operators may 
be as powerful and sensitive a harbinger of new physics as were the parton 
distributions in the mid-70 's. 

^^See also related work in Refs. |122| and |119l 112(11 11211 and for a recent review of k± 
factorization in heavy quark production, see Ref. | 12.-{| . 
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4.3 Gluon and quark production in the dense/ A A region 

{pAi/kl = pA2/kl ~ 1). 

Since pi,2/^l ~ one has thus far not been able to compute particle pro- 
duction analytically in the CGC. The problem however is well defined in 
weak coupling and can be solved numerically jl28l 11291 ll3Uj . In practice, 
this is equivalent to solving the Yang-Mills equations after the collision of 
two ultrarelativistic nuclei with initial conditions given by the classical gluon 
fields of each of the nuclei before the collision |77j . 




Figure 20: Diagrams which break k± factorization are 0(1) for kj_ < Qs-they 
are not suppressed. 

The numerical simulations performed thus far assume Gaussian initial 
conditions as in the MV model. As we noted in our discussion of Deuteron- 
Gold collisions, these may actually be good initial conditions for central 
Gold-Gold collisions at RHIC where the typical x is of order 10~^. These 
will not be good initial conditions at the LHC where the typical x at central 
rapidities will be at least an order of magnitude lower. In that case, one has 
to use solutions of JIMWLK. 

Unlike gluon production in the pp and pA cases, kj_- factorization breaks 
down in the AA-case \12S\. \127^ . k± factorization breaking diagrams of the 
sort shown in Fig. \2iA are of order 1 and there are a large number of these. 
A significant consequence is that one cannot factor the quantum evolution 
of the initial wavefunctions into unintegrated gluon distributions unlike the 
pA case. 

Nevertheless, there is a systematic way to include small x effects in the 
AA case. Numerical solutions of the JIMWLK equations now exist and one 
can study the rapidity evolution of the Wilson lines we discussed in Lecture 
II jl8j . The numerical lattice formalism developed in Ref. jl28j is ideally 
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suited to compute particle production in the forward light cone by matching 
the Wilson lines from each of the nuclei on the light cone. This program has 
not yet been carried out. 

In the following, we will restrict ourself to discussing numerical solutions 
with Gaussian initial conditions. The saturation scale Qg (which is an input 
in the numerical solutions in this approximation) and the nuclear radius R 
are the only parameters in the problem The number and energy of gluons 
released in a heavy ion collision of identical nuclei can therefore be simply 
expressed as 



1 dE 



■kB? drj g 
1 dN CN 



_2 ' 

Q's , (41) 



2 



vri?^ drj g 

where (up to 10% statistical uncertainity) ce = 0.25 and cn = 0.3. Here rj is 
the space-time rapidity. The numerical computations are performed in terms 
of the color charge squared per unit that we discussed in lecture 

II (see Eq. [TT|) . The results for ce and cat are obtained for these values. 
When re-expressed in terms of Qg, there is a small logarithmic uncertainity 
depending on the scale of Qs- 
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Figure 21: V2 from melting colored glass. 

The number distributions of gluons can also be computed in this ap- 
proach. Remarkably, one finds that a) the number distribution is infrared 
finite, and b) the distribution is well fit by a massive Bose-Einstein distri- 
bution for k±/Qs < 1.5 GeV with a "temperature" of ~ 0A7Qs and by the 
perturbative distribution Qs/k\_ for k±/Qs > 1.5. 

^"Another scale is the nucleon size-the color charge over this scale is constrained to be 
zero. Our results are insensitive to reasonable variations in this scale. 
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The RHIC data on the multiphcity (approximately 1000 hadrons in one 
unit of rapidity) and transverse energy (approximately 500 GeV for central 
rapidities) of produced hadrons combined with Eq. place strong con- 
straints on what Qs can be. If Qs is too small, we will find, absurdly, that 
the initial transverse energy is less than the final measured transverse energy. 
On the other hand, if Qs is too large, we will find that the initial multiplicity 
of gluons is greater than the final multiplicity of hadrons. Besides, the initial 
energy per particle will be too difficult to get rid of to match the oberved 
experimental value. While there is no apparent theorem that prohibits the 
initial gluon multiplicity being greater than the final hadron multiplicity, 
such a situation is unlikely in all statistical/hydrodynamic scenarios of the 
RHIC collisions. These constraints therefore allow us to place the bound 
that HSl]. 

1.3 <Qs<2 GeV (42) 

This bound is consistent with an extrapolation of the Golec-Biernat-Wusthoff 
(see Lecture I) fit of to the RHIC data (see Eq. EJ. A simple extrapola- 
tion gives Qs ~ 1.4 GeV. Other considerations from RHIC will suggest (as 
we will see shortly) that the value of Qs should be closer to the higher end 
of this bound. 
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Figure 22: The space-time evolution of a RHIC collision. 

In Ref. jl32j . we computed the elliptic flow produced by the CGC right 
after the collision. We found (see Fig. 0T|) that only about half the elliptic 
flow could be generated in this way. Besides, the p±_ dependence of the 



37 



V2{p±) distribution was much steeper than the RHIC data. This result 
strongly suggested that final state interactions were important at RHIC 
and that the RHIC data could not be explained by initial state interactions 
alone. If there are significant final state interactions and hydrodynamic 
flow (as suggested by successful hydrodynamic fits to the RHIC data), this 
suggests that Qs must be closer to the 2 GeV upper bound. Otherwise, 
if it were instead closer to the lower bound in Eq. 1411 the partonic E± at 
r] = would be uncomfortably close to the observed hadronic value. Since 
Et/N ~ 0.9 Qs, for ~ 2 GeV, the system will have to do a significant 
amount of work to reduce the transverse energy per particle to the observed 
value; this is what one expects from hydrodynamics. Thus from very simple 
considerations, and from the RHIC data, we have already learned quite 
a bit about the initial state of the matter produced at RHIC. The initial 
energy density (at a time r ~ 3/Qs) can be estimated to be between 20 — 40 
GeV/fm^, much greater than the energy density required by the lattice to 
form a Quark Gluon Plasma (QGP). 

What is not understood is the transition to the QGP from the CGC. (A 
cartoon showing the timeline of a RHIC collision is shown in Fig. 1221 ) Due 
to the rapid expansion of the system, the occupation number of modes falls 
well below one on time scales of order 1/Qs- From these times onwards, 
one expects the canonical classical approach to break down-well before ther- 
malization. On the other hand, for elliptic flow from hydrodynamics to 
be significant, the conventional wisdom is that thermalization should set 
in early. A necessary condition is that momentum distributions should be 
isotropic. The CGC initial conditions are very anisotropic with < p± >~ Qs 
and < Pz >~ 0. How does this isotropization take place? All estimates of 
final state re-scattering of partons formed from the melting CGC, both from 
2 — > 2 processes |133l 11341 1135j and 2 — > 3 processes |136j suggest thermal- 
ization takes longer than what the RHIC collisions seem to suggest-in the 
latter case, Tthermai ~ 13/5 77- , which at RHIC energies gives rthermai ~ 2 — 3 

fm [TTrTj . 

Recently, it has been suggested that collective instabilities, the non- 
Abelian analog of the well known Weibel instabilities in plasma physics, 
can speed up themalization (I38i 11391 ll4Uj . Starting from very anisotropic 
(CGC-like) initial conditions, these instabilities drive the system to isotropy 
on very short time scales, of order 1/Qs in some estimates. 

What is the relation of this language of instabilities and that of our 
classical field simulations? Why didn't we see such instabilities? One pos- 
sibility is that our particular initial conditions, the non-linearities of the 
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fields and the rapid expansion of tfie system kill the growth of instabilities. 
Another intriguing possibility is the following. In our solutions of the Yang- 
Mills equations, it is assumed that the incoming nuclei are delta-function 
sources on the light cone. This assumption directly leads to an explicit 
boost invariance of the classical fields produced in the collision. If one re- 
laxes the delta function criterion, small violations of boost invariance should 
exist. Are these the seeds for the instabilities? Further, to properly study 
thermalization, one should better understand the interaction of high mo- 
mentum (particle) and low momentum (field) degrees of freedom and their 
evolution. This leads to a real time renormalization group description jl41j 
which should be coupled with the stability analysis. Work in these directions 
is ongoing and will hopefully soon provide further insight into the nature of 
thermalization in heavy ion collisions. 

Much of the discussion about equilibration has focused on kinetic equi- 
libration. However, equally interesting is the problem of chemical equilibra- 
tion. At high energies, the initial state in a heavy ion collision is dominated 
by gluons. Sea quark pairs are produced from the gluon fields. Are they 
produced in sufficient numbers and do they re-interact sufficiently strongly 
for the system to reach chemical equilibrium (where the ratio of gluons to 
quarks is expected to be 2>2/21N f)1 Under normal circumstances, one would 
expect, in weak coupling, that the production of quarks to be suppressed. 
However, since the fields from the CGC are of order 1/g, it is conceivable 
that the strong fields could drive the system to equilibrium. First steps 
have been taken to study this problem |142| I143j . The problem involves 
solving the Dirac equation in the background field of the two nuclei (which, 
as we discussed previously, can be computed numerically). Hopefully further 
progress can be made on this problem in the near future. 

Thus far, we have said very little about heavy ion phenomenology at 
RHIC. There is a successful CGC model- the KLN model |144j -which does 
a very good job of explaining certain bulk features of the data such as the 
centrality dependence and energy dependence of the multiplicity and the 
rapidity distributions at different RHIC energies. The model assumes k±- 
factorization and parton-hadron duality. With regard to the former, as we 
have discussed, we expect A;_|_-factorization to be broken in AA-collisions. 
Thus the fact that KLN does well suggests perhaps that the effect is a 
small one, or that several factors contribute to the successful predictions. 
For instance, some effect of the factorization breaking might be taken into 
account by the 'gluon liberation' factor (our cn in Eg. I41() . What KLN has 
trouble with is the energy dependence-on face value one gets E±/N ~ Qs, 
about a factor of 3 larger than the RHIC data. This requires that post 
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CGC either the system does a lot of isentropic work or partons multiply 
furiously or some combination of the two. Further, the RHIC data on elliptic 
flow (as we have discussed) and jet quenching (which we have not) suggest 
that the produced partons re-interact strongly after the CGC stage. The 
latter was forcefully established by the RHIC Deuteron-Gold data. A model 
which assumes rapid thermalization after the CGC and subsequent ideal 
hydrodynamic evolution with KLN initial conditions was proposed by Hirano 
and Nara (145j . This apparently takes care of the E± and elliptic flow issues. 
However, several loose ends persist. For instance, even if thermalization is 
rapid, one would expect substantial entropy generation, which is not taken 
into account in matching the two descriptions. 

In hindsight, despite gaps in our understanding, the evolution in our 
understanding of heavy ion collisions has exceeded even the most optimistic 
estimates. For instance, much evidence suggests that a strongly interacting 
quark gluon fluid is briefly created at RHIC even if its properties remain 
to he fully understood The RHIC experimentalists deserve full credit 

for this by devising ingenious observables that have stretched even the most 
fortunate of models, with most falling like nine pins. Experiments are driving 
theory, disciplining our fancy, and that's how it should be. 
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